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Abstract 

The free energy for QCD at high temperature T is calculated to order g 5 using 
effective-field-theory methods to separate the contributions from the momentum scales 
T and gT. The effects of the scale T enter through the coefficients in the effective 
lagrangian for the 3-dimensional effective theory obtained by dimensional reduction. 
The perturbation series for these coefficients seem to be well-behaved if the running 
coupling constant is sufficiently small: a s (2irT) -C 1. For the contribution to the free 
energy from the scale gT, the perturbation series is well-behaved only if a s (2irT) is an 
order of magnitude smaller. The implication for applications of perturbative QCD to 
the quark-gluon plasma are briefly discussed. 
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One of the most dramatic predictions of quantum chromodynamics (QCD) is that when 
hadronic matter is raised to a sufficiently high temperature or density, it will undergo a 
phase transition to a quark-gluon plasma. One of the major thrusts of nuclear physics 
in the next decade will be the effort to study the quark-gluon plasma through relativistic 
heavy-ion collisions. For this effort to be successful, it will be important to understand 
the properties of the plasma as accurately as possible. The two major theoretical tools for 
studying the quark-gluon plasma are lattice gauge theory and perturbative QCD. Lattice 
gauge theory has the advantage that it is a nonperturbative method and applies equally 
well to the hadronic phase. It is an effective method for calculating the static equilibrium 
properties of a plasma with baryon density, but it can not be easily applied to problems 
involving dynamical properties or to a plasma that is away from equilibrium or has nonzero 
baryon density. These problems can however be studied using perturbative QCD, provided 
that the temperature T of the plasma is sufficiently high. As the temperature decreases, the 
running coupling constant g of QCD increases, causing perturbation theory to break down 
at some temperature above the critical temperature T c for the phase transition. One of the 
basic questions in the theory of the quark-gluon plasma is how large must T be in order 
for perturbative QCD to be applicable. Is this method useful at temperatures that may be 
achievable in heavy-ion collisions, which are at most several times T c l 

In order to answer this question, it is necessary to understand the structure of the pertur- 
bation series to all orders and also to carry out explicit higher-order calculations. The first 
step has been carried out for the free energy ffl. The structure of the perturbation series 
is nontrivial, because a strict perturbation expansion in g 2 has severe infrared divergences 
associated with the exchange of static gluons. Physically, these divergences are screened 
by plasma effects. The screening of electrostatic gluons can be taken into account by a 
resummation of perturbation theory, but the screening of magnetostatic gluons can only be 
treated using nonperturbative methods. Once the structure of the perturbation expansion is 
understood, it is still necessary to carry out explicit perturbative calculations to determine 
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quantitatively how high the temperature must be in order for perturbative QCD to be ac- 
curate. Only recently has the calculational technology of thermal field theory progressed to 
the point where it is possible to carry out perturbative calculations to a high enough order 
that the running of the coupling constant comes into play [2-8]. The first such calculation 
for the quark-gluon plasma was the calculation of the free energy to order g 4 by Arnold and 
Zhai in 1994 M. This calculation has recently been extended to order g 5 by Kastening and 
Zhai U and by the authors ||. In this Letter, we summarize the calculations of Ref. [H] and 
discuss implications for the application of perturbative QCD to the quark-gluon plasma. 

The static equilibrium properties of a quark-gluon plasma at temperature T are governed 
by the free energy density F = — (T/V) log^QCD, where V is the volume of space. The 
partition function Zqcd is given by a functional integral over quark and gluon fields on a 
4-dimensional Euclidean space-time, with the Euclidean time r taking its values on a circle 
with circumference 1/T. In the limit in which the quarks are massless, the free energy is a 
function of T and a s = g 2 /(47r) only. 

The structure of the weak-coupling expansion for the free energy to all orders was deduced 
in Ref. PJ. The free energy can be expressed as the sum of three contributions coming from 
the momentum scales T, gT, and g 2 T: 

F = [f E (A E ) + /m(Ae,A m ) + f G (A M )]T, (1) 

where is an arbitrary factorization scale that separates the scales T and gT, while Am 
separates the scales gT and g 2 T. The contributions from the three momentum scales can 
be unraveled by constructing a sequence of two effective field theories. The first effective 
theory, electrostatic QCD (EQCD), is a 3-dimensional Euclidean field theory involving the 
electrostatic gauge field Ao( x ) an d the magnetostatic gauge field A*(x). The lagrangian for 
EQCD is 

£eqcd = \G%G% + ^(AA ) a (AA ) a + \m\AlAl + 5£ EQCD , (2) 
where Di is the covariant derivative for the adjoint representation with coupling constant g E 
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and G°j is the magnetostatic field strength. The term <5£eqcd in (@) includes all other local 
gauge-invariant operators that can be constructed out of A and Ai, including nonrenor- 
malizable interactions. Static gauge-invariant correlation functions in thermal QCD can be 
reproduced at long distances by tuning g E , m E , and the parameters in <5£eqcd as functions 
of g, T, and the ultraviolet cutoff A E of EQCD. In physical quantities, the A B -dependence 
of the parameters cancels the A^-dependence from loop integrals in EQCD. Since the pa- 
rameters of EQCD take into account effects from the scale T, they can be calculated as 
perturbation series in the running coupling constant a s (fi) with \x of order T. 

Since the free energy is a static quantity, it can be calculated using EQCD. The free 
energy can be written 



where ^eqcd is the partition function for EQCD and f E is the coefficient of the unit operator, 
which was omitted from the effective lagrangian (|2|). This coefficient gives the contribution 
to the free energy ([!]) from the momentum scale T. The logarithm of -Zeqcd i n © includes 
the contributions fu and fa from the scales gT and g 2 T, respectively. These contribu- 
tions can be separated by constructing a second effective field theory, magnetostatic QCD 
(MQCD), which involves only the magnetostatic gauge field A"(x). The term ju in ([!]) is 
the coefficient of the unit operator in the lagrangian for MQCD. It can be computed using 
perturbative methods as an expansion in powers of g starting at order g 3 . The term fa 
in ([!]) is proportional to the logarithm of the partition function of MQCD. It can only be 
calculated using nonperturbative methods. Surprisingly, it can be expanded in powers of g 
beginning at order g 6 , with coefficients that can be calculated using lattice simulations of 
MQCD 0. Since we only calculate the free energy to order g 5 , we do not consider the term 
fa any further. 

To calculate the free energy to order g 5 , the only parameters of EQCD that are required 
are g 2 E to leading order in g 2 and m 2 E and f E to order g A . The gauge coupling constant g E for 
EQCD is determined at leading order simply by comparing the lagrangians for EQCD and full 




(3) 
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QCD: g\ = g 2 T. The other two parameters can be determined by computing static quantities 
in both full QCD and EQCD, and demanding that they match. It is convenient to carry 
out these matching calculations using a strict perturbation expansion in g 2 . This expansion 
is afflicted with infrared divergences due to long-range forces mediated by static gluons, 
and an infrared cutoff is therefore required. Physically, these divergences are screened by 
plasma effects, but screening is not taken into account in the strict perturbation expansion. 
Nevertheless, this expansion can be used as a device for determining the parameters of 
EQCD, since they depend only on short distances of order 1/T. 

The parameter rriE can be determined by matching the strict perturbation expansions 
for the electric screening mass m e \ in full QCD and in EQCD. Beyond leading order in g, m e \ 
becomes sensitive to magnetostatic screening and requires a nonperturbative definition [0. 
However, in the presence of an infrared cutoff, the electric screening mass can be defined in 
perturbation theory by the location of the pole in the propagator for A (r, x) at spacelike 
momentum (ko = 0, k). Denoting the appropriate component of the gluon self-energy tensor 
by U(k 2 ) where k 2 = k 2 , we must solve the equation k 2 + U(k 2 ) = at k 2 = — to^. 
Since the solution rn? eX is of order g 2 , we can expand U(k 2 ) as a Taylor series around k 2 = 0. 
To determine m\ to order g 4 , we must calculate 11(0) to two-loop accuracy and n'(0) to 
one-loop accuracy. We use dimensional regularization with 3 — 2e spacial dimensions to 
cut off both infrared and ultraviolet divergences. The sums and integrals can be evaluated 
analytically using methods developed by Arnold and Zhai M. The resulting expression for 
is an expansion in integral powers of a s . There is no aPj 2 term, unlike in the expression 
for that correctly incorporates the effects of electrostatic screening ||. This g 3 term 
arises because the g 4 correction includes a linear infrared divergence that is cut off at the 
scale gT. Since we use dimensional regularization as an infrared cutoff, this power infrared 
divergence is set equal to 0. 

In EQCD with an infrared cutoff, the electric screening mass m e \ can be defined in 
perturbation theory by the location of the pole in the propagator for the field A (x). 
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Denoting the self-energy function for A (x) by ITe(A; 2 ), the screening mass m e \ satisfies 
k 2 + m 2 E + IIe(A; 2 ) = at k 2 = — m 2 y In the strict perturbation expansion for EQCD, we 
treat perturbation parameter of order g 2 . After Taylor-expanding ITe(A; 2 ) around 

k 2 = 0, there is no scale in the loop integrals, so they all vanish with dimensional regular- 
ization. The solution for the screening mass is therefore trivial: m 2 j = m 2 E . Matching this 
result with the strict perturbation expansion from the full theory and taking the limit e — > 0, 
we find 
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e=0 



4vra s (/i)T 2 |l + |n F + 0.612 - 0A88n F - 0.0428 



2nT 
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where rip is the number of flavors of quarks and /i is the renormalization scale for the QCD 
coupling constant. The order-e terms in m 2 E are also required in the calculation of the free 
energy. These terms are given by 



dm 2 E 



= y 

e=0 



,W3.97 + 21ogA| + (o.597 +i log ^)nA (5) 



where the infrared cutoff Ap is the scale introduced by dimensional regularization. 

The coefficient can be determined by matching the strict perturbation expansions for 
the free energy in full QCD and in EQCD. In full QCD, the free energy F is calculated to order 
g 4 by evaluating the sum of vacuum diagrams through three-loop order, using dimensional 
regularization to cut off both infrared and ultraviolet divergences. The resulting expression 
for F is an expansion in integral powers of a s . There is no o^J 2 term, in contrast to the 
expression for the free energy that correctly includes the effects of electrostatic screening 
0, ID]. This g 3 term arises because the g A correction includes a linear infrared divergence 



that is cut off at the scale gT. In the strict perturbation expansion, this term appears as a 
power infrared divergence that is set to zero in dimensional regularization. 

In EQCD, the free energy is given by (^). All the loop diagrams in the strict perturbation 
expansion for log^EQCD vanish with dimensional regularization, since there is no scale for 
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the integrals. The only contribution to logi^EQCD comes from the counterterm SJ'e which 
cancels logarithmic ultraviolet divergences proportional to the unit operator. The resulting 
expression for the free energy is simply F = (/# + 5/e)T. The counterterm is determined 
by calculating the ultraviolet divergent terms in logZ E Q CD . If we use dimensional regular- 
ization together with a minimal subtraction renormalization scheme in EQCD, then SJe is a 
polynomial in g E , m 2 E , and the other parameters in the EQCD lagrangian. The leading term 
in 8/e is proportional to g E m E , and its coefficient can be determined by a simple 2-loop 
calculation: 

3 

6 fx = —^9 2 E m E- (6) 

When this counterterm is expressed in terms of the parameters g and T of the full theory, 
we must take into account the fact that m 2 E in multiplies a pole in e. Thus, in addition 
to the expression for m E given in (|j), we must also include the terms of order e which are 
given by (^). Matching F = (/# + 8/e)T with the strict perturbation expansion for F in 
the full theory, we obtain 



/,(A,) = -^{l + l»,-H( 1 + A n ,)^0 + 

U,;j fl + s^) ( x - 1^) lo s Sf - 135 ( x + loF 



244.9 - 17.24n F - 0.415n| 



a - 2 



2nT — v - 1 e -J-~°2tiT\ Vtt ) J" 

This expression differs from that given in Ref. |l|], where the counterterm ([|) was not taken 
into account. 

We have calculated two terms in the perturbation series for m 2 E and three terms in the 
series for We can use these results to study the convergence of perturbation theory for 
the parameters of EQCD. We consider the case of he = 3 flavors of quarks, although our 
conclusions will not depend sensitively onnp. The question of the convergence is complicated 
by the presence of the arbitrary renormalization and factorization scales \i and A#. The next- 
to- leading-order (NLO) correction to Je is independent of \i and A E , and is small compared to 
the leading-order (LO) term provided that a s (/i) < 1.1. The NLO correction to m E and the 
next-to-next-to-leading-order (NNLO) correction to Je both depend on the renormalization 
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scale jj,. One scale-setting scheme that is physically well-motivated is the BLM prescription 



II], in which \i is adjusted to cancel the highest power of rip in the correction term. This 
prescription gives fi = 0.93 nT when applied to m E and /i = 4.4 nT when applied to 
These values differ only by about a factor of 2 from 2nT, which is the lowest Matsubara 
frequency for gluons. Below, we will consider the three values fi = nT, 2nT, and 4nT. For 
the NLO correction to m E to be much smaller than the LO term, we must have ct s (fi) -C 0.8, 
3.8, and 1.4 if /x = nT, 2tiT, and 4nT, respectively. Based on these results, we conclude 
that the perturbation series for the parameters of EQCD are well-behaved provided that 
a s (2nT) < 1. 

The NNLO correction for depends not only on fi, but also on the factorization scale 
Ae- Because the coefficient of log(A£i/27rT) in (|^) is so much larger than that of \og(fi/2nT), 
the NNLO correction for is much more sensitive to Ae than to fi. It is useful intuitively to 
think of the infrared cutoff as being much smaller than the ultraviolet cutoff fi. However, 
these scales can be identified with momentum cutoffs only up to multiplicative constants 
that may be different for fi and Ae- Both parameters are introduced through dimensional 
regularization, but /i arises from ultraviolet divergences of 4-dimensional integrals, while 
Ae arises from infrared divergences of 3-dimensional integrals. We might be tempted to set 
Ae = fi, but then the NNLO coefficient in is large. For the choice fi = 2ttT, the correction 
to the LO term is a multiplicative factor 1 — 0.9a s + 6AQa 2 . The NNLO correction can be 
made small by adjusting A#. It vanishes for A# = 5.8 tcT, 5.1 ttT, and 4.5 7rT if fi = nT, 
27rT, and 4irT, respectively. We conclude that the perturbation series for Je is well-behaved 
if the factorization scale A# is chosen to be approximately 5nT. Whether this choice is 
reasonable can only be determined by calculating other EQCD parameters to higher order 
to see if the same choice leads to well-behaved perturbation series. 

The choice of A^ that makes the perturbation series for the EQCD parameters well- 
behaved may be much larger than the largest mass scale m# of EQCD. Perturbative correc- 
tions in EQCD will then include large logarithms of Ae/tue- This problem can be avoided 
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by using renormalization group equations to evolve the parameters of EQCD from the initial 
scale A down to some scale A' E of order m#. The coefficient satisfies the renormalization 
group equation 

A B — — f E = ^—q%m%. (8) 

dA E 2tt 2HE e v ' 

The evolution of g E and m E occurs only at higher order in the coupling constant and therefore 
can be ignored. The solution to the renormalization group equation is therefore trivial: 

MA*) = f E (A E ) - T^glml log ^ . (9) 

Having determined the parameters of EQCD to the necessary accuracy, we proceed to 
calculate the free energy using (fj). The contribution from the scale T is given by the 
coefficient in (0). The contribution from the scale gT is given by fu = — log^EQCD/^- 
In order to calculate ju using perturbation theory in EQCD, we must include the effects of 
the mass parameter m 2 E to all orders, but the gauge coupling constant gE can be treated as 
a perturbation parameter. The contributions to log^EQCD of orders g 3 , g A , and g 5 are given 
by the 1-loop, 2-loop, and 3-loop vacuum diagrams in EQCD, respectively. The integrals 



can be calculated analytically using methods developed by Broadhurst |12| . The two- loop 
integrals include an ultraviolet pole in e that is proportional to g E m E . This divergence is 
cancelled by the counterterm SJe for the coefficient of the unit operator, which is given in 
(|). Our final result for the coefficient fu in (|l|) is 

^2 / „2 \ 2 



f M (A E ) = ~izm\ 



2 

3?r' 



0.256 -- log A^)^L_ _ 27.6 (-**- 

4 vriE/ 2nrriE \2iimi 



(10) 



Note that the dependence of /a/ on A £ cancels that of in (|^). The expression ( |T0D can 
be expanded in powers of g by setting g E = g 2 T and using the expansion ([|) for m E . 

We now consider the convergence of the perturbation series ( |10D for Jm- The size of 
the NLO correction depends on the choice of the factorization scale A#. It is small if Ae 
is chosen to be approximately uie- The NNLO correction in (|T0D is independent of any 
arbitrary scales. If rip — 3, it is small compared to the leading order term only if a s -C 0.17. 
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Thus the perturbation series for Jm is well-behaved only for values of a s (2irT) that are much 
smaller than those required for the parameters of EQCD to have well-behaved perturbation 
series. 

Adding (|7]) and (10) and expanding in powers of v / al, the complete expression for the 
free energy F is 



8^ 
' 45 



rp4 



3/2 



7T 



+ F 4 

7T / V 7T 



Q', 



+ F 5 



a. 
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where the truncation error is of order log a s . The coefficients in this expansion are 



F 2 
F 3 
F, 



l + §n F 



15 



30 fl + \n F 



3/2 



237.2 + 15.97n F - 0.413n^ + (l + \n F \ log ( — (l + \n F 

2 V 7T 



165 



1 + &n F ) (1 - in F ) log 



33 



2nT ' 



(12) 
(13) 
(14) 

(15) 



(i + i 



Tip 



1/2 



- 799.2 - 21.96n F - 1.926n| 



+ ^ (l + |n F ) (l - ±n F ) log 



2ttT 



(16) 



The coefficient _F 2 was first given by Shuryak ||, while F 3 was first calculated correctly 
by Kapusta [JTIJ. The coefficient _F 4 was computed in 1994 by Arnold and Zhai ||. The 
coefficient F 5 in ( |T6"D has been calculated independently by Kastening and Zhai using a 
different method ||. 

We now ask how small a s must be in order for the expansion ( |TT| ) to be well-behaved. 
For simplicity, we consider the case n F = 3, although our conclusions are not sensitive 
to np. If we choose the renormalization scale \x = 2nT motivated by the BLM criterion 
H|, the correction to the LO result is a multiplicative factor 1 — 0.9a s + 3.3a^ 2 + (7.1 + 



3.51oga s )a 2 — 20. 8a F ^ 2 . The oC'J 2 term is the largest correction unless a s (27rT) < 0.12. We 
can make the a b J 2 term small only by choosing the renormalization scale to be near the value 
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fj, = 36.57rT for which F 5 vanishes. This ridiculously large of fi arises because the scale fi has 
been adjusted to cancel the large g 5 correction to ju in (P~0|) . This contribution arises from 
the momentum scale gT and has nothing to do with renormalization of a s . We conclude 
that the expansion (|TTD for F in powers of ^/a^ is well-behaved only if a s (27rT) is an order 
of magnitude smaller than the value required for the EQCD parameters to be well-behaved. 

We now consider briefly the implications for theoretical studies of the quark-gluon plasma. 
We have found that the convergence of perturbation theory requires much smaller values of 
a s {2'KT) for quantities at the scale gT than for quantities at the scale T. The critical 
temperature T c for formation of a quark-gluon plasma is approximately 200 MeV. It may be 
possible in heavy-ion collisions to produce a quark-gluon plasma with temperatures several 
times T c . At T = 350 MeV, a s (2irT) pa 0.3, which is small enough that perturbation theory 
may be reasonably convergent at the scale T, but it is certainly not convergent at the scale 
gT. We conclude that at the temperatures achievable in heavy-ion collisions, perturbative 
QCD may be accurate when applied to quantities that involve the scale T only. However 
nonperturbative methods are required to accurately calculate quantities that involve the 
scales gT and g 2 T associated with screening in the plasma. 
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